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ABSTRACT 

We propose use of the thickness of the ionisation front as a discriminant between alterna- 
tive modes of reionisation in the early universe, by stars or by miniquasars. Assuming a 
photoionisation-recombination balance, we find that for miniquasar sources the transition 
from neutral to ionised intergalactic medium is extended and has two features. The first is 
a sudden steep increase in the neutral fraction with a typical width of 5 — 10 comoving mega- 
parsecs, depending on the miniquasar power. The second feature is a long wing that represents 
a much slower transition from neutral fraction of m . 8 to 1 . The angular resolution of LOFAR 
is expected to resolve these scales and will, therefore, play an important role in discriminating 
the hard sources of ionising photons from the stellar ones. 

Key words: galaxies: cosmology: theory - large-scale structure of Universe - diffuse radia- 
tion - radio lines: general - quasars: general 
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1 INTRODUCTION 

The universe is generally considered to have become reionised at 
a redshift larger than 10. This result, coming from analysis of the 
WMAP polarisation power spectrum (Kogut et al. 2003 & Spergel 
et al. 2003), came as a surprise, in the context of earlier studies of 
the Gunn-Peterson effect inferred to be present at z ~ 6 (e.g., Fan, 
X., et al. 2002) and of the high temperature of the intergalactic 
medium at z ~ 3 (Theuns et al. 2002; Hui & Haiman 2003). 
Given this result, the ionising sources cannot be known quasars 
or normal galaxies. Rather, recourse must be had to Population III 
stars or to miniquasars, both of which represent hypothetical but 
plausible populations of the first objects in the universe, and that 
are significant sources of ionising photons. 

The impact of stellar sources have been studied by many 
authors (e.g., Cen 2003, Ciardi, Ferrara, & White 2003, Haiman 
& Holder 2003, Wyithe & Loeb 2003, Sokasian et al. 2003, 
Somerville & Livio 2003). These studies find in general that in 
order to provide enough ionising flux at or before z = 15, for 
the usual scale-invariant primordial density fluctuation power spec- 
trum, one needs Population III stars, which provide about 20 times 
more ionising photons per baryon than Population II (Schaerer 
2002, Bromm, Kudritzki, & Loeb 2001), or an IMF that initially is 
dominated by high mass stars (Daigne et al. 2004). This is in agree- 
ment with recent numerical simulations of the formation of first 
stars from primordial molecular clouds that suggest that the first 
metal-free stars were predominantly very massive, m* £ IOOA/q 
(Abel, Bryan, & Norman 2000; 2002 & Bromm, Coppi & Larson 
2002). 

Theoretically, there is some tension between the low ampli- 
tude of the fluctuations as interpreted from low redshift data (e.g., 
Croft et al. 2002, Zaroubi et al. 2004) and the early ionisation ob- 



served by the WMAP satellite. Most of the theoretical and numeri- 
cal results mentioned earlier require a somewhat higher value of the 
fluctuation amplitude and would have a harder time satisfying all of 
the observational constraints with lower amplitude (trg ~ 0.8). It 
is unclear whether this tension is caused by the theoretical uncer- 
tainties in the physical details of the ionisation process or it is an 
indication of the existence of a non-standard cosmological model 
in which nonlinear structures assemble earlier than in the normal 
ACDM scenario (e.g., Sugiyama, Zaroubi & Silk 2004; Avelino & 
Liddle 2004). 

Miniquasars have also been considered as a significant ion- 
ising source (e.g., Ricotti & Ostriker 2004a; 2004b, Madau et al. 
2004, Oh 2000, 2001, Dijkstra, Haiman & Loeb 2004). These lat- 
ter are as plausible as Population III stars, whose nucleosynthetic 
traces have not yet been seen even in the most metal-poor halo stars 
nor in the high z Lyman alpha forest, in view of the correlation 
between central black hole mass and spheroid velocity dispersion 
(Ferrarese 2002, Gebhardt et al. 2000). This correlation demon- 
strates that seed black holes must have been present before spheroid 
formation. Indeed, recent observations of a quasar host galaxy at 
z = 6.42 (Walter et al. 2004) (and other quasars) suggest that su- 
permassive black holes were in place and predated the formation of 
the spheroid. Theory suggests that the seeds from which the SMBH 
formed amounted to at least IOOOM0 and were in place before 
z ~ 10 (Islam, Taylor & Silk 2003, Madau & Rees 2001). 

Can we distinguish between the alternative hypotheses of stel- 
lar versus miniquasar ionisation sources? One distinguishing fea- 
ture is the intrinsic source spectrum, which is thermal for stars but 
with a cut-off at a few times the Lyman limit frequency, but power- 
law for miniquasars with a spectrum that extends to high energies 
with nearly equal logarithmic increments in energy per decade of 
frequency. We show here that there is a dramatic difference between 
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these two cases in the widths of the ionisation fronts. Only the mini- 
quasar model translates to scale-dependent 21 cm brightness tem- 
perature fluctuations that should be measurable by LOFAR. 



2 WIDTH OF X-RAY IONISATION FRONTS 

Typically, the spectrum of the UV radiation of stellar sources is 
roughly flat with a cutoff at a few times the Lyman limit frequency. 
This results in a very sharp transition from ionised to neutral IGM. 
The thickness of such a region may be estimated from the mean 
free path of a photon with energy E within a neutral medium, 



n H OH{E) ' 



(1) 



where n H ~ 2.2 x 10" 7 cm~ 3 (l + z) 3 (Bennett et al. 2003) 
is the mean number density of hydrogen at a given redshift, and 
<th(E) = (To (Eo/E) is the bound-free absorption cross-section 
for hydrogen with a = 6 x 10~ 18 cm 2 and E = 13.6 eV . At 
z = 10, the mean free path for the most energetic photon from a 
stellar source is of the order of 0.1 comoving kpc. 

In contrast, the power-law behavior of the radiation specific 
intensity of miniquasars produces a much thicker ionisation front. 
In the following calculation, we estimate this thickness assuming 
a uniform mass density distribution in which a miniquasar is em- 
bedded. This obviously does not take into account any mass or 
quasar clustering properties; however, for the purpose of explor- 
ing the thickness to within an order of magnitude and its qualitative 
evolution with redshift and dependence on the black hole mass, this 
assumption is acceptable. A more realistic numerical calculation 
will be carried out in the future. 

We are also going to assume ionisation-recombination equilib- 
rium. This may be safely assumed for the photons we are interested 
in as they have a mean free path of ~ 1 — 10 comoving Mpc at red- 
shift ~ 10 — 5, which translates to a mean scattering time of about 
0.5-5 million years, much less than a Hubble time at the relevant 
redshift range. 

The miniquasar source spectrum is taken to be 

7V en A TP — 0.8 -E/E c -2 -1 

1(E) = AgE e ' cm s 



(2) 



where g is the Gaunt factor and A is normalized such that the lumi- 
nosity of the miniquasar is 10 43 ^i erg s _1 , for black holes of mass 
10 5 /uM©, with 1 < < 10 3 , and the cutoff energy E c is 200 keV 
(Sazonov et al. 2004). This translates to a number of ionisations at 
a distance r from the source, 



N(E;r) = e - T(E '- r) 
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with 



r(E;r)= I n H x H io(E)dr. (4) 
Jo 

Here xhi is the hydrogen neutral fraction. 

The factor E/Ei in equation [3] is added to account for the 
mean number of ionisations per ionising photon with Ei ~ 36 eV. 
This last choice is a simplifying approximation to the number of 
ionisations per photon which has a more complex dependence on 
energy and ionised fraction of hydrogen (Shull & van Steenberg 
1985; Dijkstra, Haiman & Loeb, 2004). The approximation we 
adopt here is reasonable for photons within the energy range that 
gives rise to the thickness of the ionisation front, i.e., photons with 
E « 0.2 - 1 keV. 



The cross-section quoted earlier does not take into account 
the presence of helium. In order to include the effect of helium, we 
follow Silk et al. (1972) who modified the cross section to become 



ct(E) = <th(E) 



n-He 



(f) 



(5) 



A proper treatment of the effect of helium is accounted for by defin- 
ing <ti to be a step function at the two helium ionisation energies 
corresponding to He I and He II. This however includes lengthy cal- 
culations and complicates the treatment, and we therefore choose 
(Ji to be a smooth function of E, an approximation that will overes- 
timate <j(E) for low energy photons. For photons with the relevant 
energy range, namely those responsible for creating the thick tran- 
sition layer, this is a reasonable assumption. 

The equation of ionisation equilibrium which expresses the 
balance between the rate of number density of ionisation and re- 
combination events is, 



a^%(l - xhi) 2 = T(E; r) n H x H i ( 1 + 



<JH tlH 



(6) 



Here Y(E\ r) is the ionisation rate per hydrogen atom for a given 
photon energy at distance r from the source. Since we are interested 
in the detailed structure of the ionisation front, Y is calculated sep- 
arately for each r using the expression, 



T(E-r) = 
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where we have introduced the following simplifying notation: 
equation |4] is written as r(E;r) = ro(r) (Eo/E) 3 where ro(r) 
is the optical depth for photons with the Lyman limit frequency at a 
given r. We also changed the integration variable to y = (Eq /E) 3 . 



.,(2) 



is the recombination cross-section to the second excited 



atomic level and has the values of 2.6 x 10 13 T 4 0,85 cm 3 s \ 
with T4, the gas temperature in units of lOOOO-R', assumed here 
to be unity. 

FigureQshows the solution of equation|6|for miniquasars with 
various luminosities and different redshifts. The transition from 
neutral to ionised IGM is extended and has two features. The first 
is a sudden steep increase in the neutral fraction with a typical 
width of 5 — 10 comoving megaparsecs (see left panel). The sec- 
ond feature is a wing that represents a much slower transition from 
xhi ~ 0.8 to 1. In the next section, we will show that with LOFAR, 
one should be able to resolve both features for the luminous mini- 
quasars. For the less luminous quasars, one might have a chance of 
just detecting the wing feature with LOFAR. Only a radio telescope 
like SKA (the Square Kilometer Array) that has 10 times the total 
collecting area of LOFAR will be able to see the whole transition 
down to miniquasars with 10 43 erg s _1 luminosity. 



3 LOFAR PREDICTIONS 

In radio astronomy, where the Rayleigh-Jeans law is usually ap- 
plicable, the radiation intensity, I(y) is expressed in terms of the 
brightness temperature, so that 



2v 2 

i(u) = ^-m 



(10) 
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Figure 1. The left panel shows the H I fraction as a function of comoving distance from the source for redshift values of z = [7.5, 10, 12.5, 15] for a source 
with fi = 10 3 in units of the Eddington luminosity for a 10 5 Mq black hole. The right panel is similar to the left one except that here the redshift is fixed to 
z=7.5 and the various curves reflect values of fi = 1, 10, 100, 1000. 



where v is the radiation frequency, c is the speed of light and k is 
Boltzmann's constant (Rybicki & Lightman 1979). This in turn can 
only be detected differentially as a deviation from Tqm b , the cos- 
mic microwave background temperature. The predicted differential 
brightness temperature deviation from the cosmic microwave back- 
ground radiation is given by (Field 1958; 1959; Ciardi & Madau 
2003), 




Here T s is the spin temperature and O m and flf, are the mass and 
baryon density in units of critical density. In the following calcula- 
tion, we take T a to be significantly larger than Tcm b ■ The enhance- 
ment in T a relative to Tqmb could be caused by the X-ray photons 
of the miniquasar itself or by an X-ray background produced col- 
lectively by miniquasars formed at higher redshifts (Ricotti & Os- 
triker, 2004a; Nusser 2004, Tozzi et ah 2000). Furthermore, we 
are adopting a standard universe with a flat geometry, Qt,h = 0.02, 
Q m = 0.3 and fl A = 0.7. 

Currently, there are a few experiments (e.g., LOFAR, PAST 
and MWA) that are being designed to directly measure 8Tt of the 
HI 21 cm hyperfine line and probe the physics of the reionisation 
process by observing the neutral fraction of the IGM as a function 
of redshift. In this study, we focus on predictions for LOFAR, but 
our conclusions could be easily applied to the other telescopes. 

The LOFAR array consists geographically of a compact core 
area and 45 remote stations. Each remote station will be equipped 
with 100 High Band antennas, 100 Low Band antennas. In the core 
area, with 2 km diameter, there will probably be 32 substations. 
For the astronomy application, there will be a total of 3200 High 
Band and 3200 Low Band antennas in the core area. Currently, the 



planned maximum baseline between the remote stations is roughly 
100 km. The Low Band antenna will be optimised for the 30-80 
MHz range while the High Band antenna will be optimised for the 
120-240 MHz range. 

The High Band antennas are sensitive enough to allow for the 
detection of the brightness temperature produced by the high red- 
shift 21 cm transitions. In this band, LOFAR will have the ability to 
measure brightness fluctuations as low as 5 mK with spectral res- 
'^iiution of ~ 1 MHz and spatial resolution of about 3 arcmin. As 
currently configured, LOFAR will be sensitive to the 21 cm emis- 
sion at redshift 6-1 1.5 over a field of 100 sq degrees. First operation 
is foreseen towards the end of 2006. 

The resolution of LOFAR depends on the telescope baseline 
and observed wavelength. As currently designed, the epoch of ion- 
isation signal will be sensitively measured mainly in the LOFAR 
core which has a FWHM resolution of ~ 3 arcmin. For more infor- 
mation, see the LOFAR web site: www.lofar.org 

Here we have convolved the predicted brightness temperature 
fluctuations with the LOFAR beam, taking a Gaussian beam of 
FWHM 3 arcmin, and obtained the results shown in figure|2| This 
figure shows the expected LOFAR observed differential brightness 
temperature for various values of redshifts with fixed central black 
hole mass of lO 8 A/0 (left panel) and for different values of black 
hole masses with fixed redshift (right panel). It is clear that the 
ionisation front for the fi = 100 — 1000 case is larger than the 
beam FWHM and could be easily observed with LOFAR. How- 
ever, this is not the case for a miniquasar with central black hole 
mass ^ 10 6 Mq. The LOFAR baseline will allow higher resolu- 
tion observations but for those baselines the collecting area of the 
telescope with respect to the core increases very little, e.g., a base- 
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Figure 2. The left panel shows 57), convolved with the 3 arcmin FWHM LOFAR Gaussian beam as a function of angle from the source center in arcmin for 
redshift values of z = [7.5, 10, 12.5, 15] for a source with [i = 10 3 corresponding to a central black hole mass of 10 8 Mq. The right panels are similar 
to the left one except that here the redshift is fixed to z=7.5 and the various curves reflect values of fi = 1, 10, 100, 1000. The ionisation front convolution 
corresponding to a Stromgren sphere, e.g., sharp ionisation front appropriate for (stellar) thermal sources, at z=7.5 is shown as well. 



line of 5 km has about 1.2 the core collecting area, resulting in a 
sharp decrease in the telescope sensitivity at these resolutions. To 
fully resolve such ionisation fronts, one would have to wait, e.g., 
for the Square Kilometer Array which for the same baseline has an 
order of magnitude higher sensitivity than LOFAR. The extended 
wings of the ionised region shown in figures □ & |2| could also 
be observed with LOFAR, possibly even for the low power mini- 
quasars. These correspond to a small but not negligible ionised frac- 
tion (xhi ~ 0.8) with a 8Tb of the order of a few. The extended re- 
gion of these low ionisation wings compensates for their low signal 
and might possibly render them observable with LOFAR. 



4 SUMMARY 

In this paper we have shown that for very simple assumptions, LO- 
FAR could in principle indicate whether a certain ionisation front 
is caused by miniquasars or stars by measuring its width. Our as- 
sumptions include a single spherical symmetric ionisation source 
with no overlapping ionised spheres, a single ionising population 
and a homogeneous background with no clustering. In the future, 
we will include the source clustering effects and the effects of the 
source luminosity function, and compute the LOFAR correlation 
function. 

In reality, however, the reionisation process might have been 
much more complex and have been induced by multiple types of 
ionisation radiation sources that were dominant at various stages of 
the universe history and present in an inhomogeneous medium. For 
example, massive Population III stars might die early on and cre- 
ate miniquasars that can maintain the IGM to be partially ionised 



through their hard photons (Ricotti & Ostriker 2004a; 2004b; 
Madau et al. 2004) while the full ionisation could be attributed 
to the soft miniquasar photons at later stages. Notwithstanding the 
complex ionisation history of the IGM, the thickness of the ionisa- 
tion fronts expected to be observed via the 21 cm experiments e.g., 
LOFAR and PAST (Peterson, Pen and Wu, 2004) should still be a 
very strong discriminator between stellar and accreting black hole 
ionising radiation sources. 
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